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We derive the general form of a master equation describing the interaction of an arbitrary multipartite quan- 
tum system, consisting of a set of subsystems, with an environment, consisting of a large number of sub- 
envirobments. Each subsystem "collides" with the same sequence of sub-environments which, in between the 
collisions, evolve according to a map that mimics relaxations effects. No assumption is made on the specific 
nature of neither the system nor the environment. In the weak coupling regime, we show that the collisional 
model produces a correlated Markovian evolution for the joint density matrix of the multipartite system. The 
associated Linblad super-operator contains pairwise terms describing cross correlation between the different 
subsystems. 
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In the study of the open dynamics of a multipartite quantum 
system S the simplifying assumption that each subsystem in- 
teracts with its own local environment is frequently made. In 
quantum communication 1 1 1 where S is identified with the set 
of information carriers employed in the signaling process, this 
is equivalent to saying that a given communication channel 
is memoryless i.e. that it acts independently on each sepa- 
rate carrier. In recent years, however, the study of correlated 
channels - sometimes called also channels with memory - has 
shown that interesting new features emerge when one makes 
the realistic assumption that the action of the noise tampering 
with the communication line is correlated over consecutive 
carriers (e.g. see f2-8| and references therein). Such corre- 
lations have been phenomenologically described in terms of 
a Markov chain which gives the joint probability distribution 
of the local Kraus operators acting on the elements of 5 |2|. 
Alternatively they have been effectively represented in terms 
of local interactions of the carriers with a common multipar- 
tite environment which is originally prepared into a correlated 
(possibly entangled) initial state ||6|, or with a structured envi- 
ronment composed by local and global components 1*3^-51. 

The aim of the present paper is to provide a continuous 
time description of correlated quantum channels in terms of 
a joint Master Equation (ME) ||9l [lOl for S. This will lead 
us to identify the structure of the Lindblad generators which 
are responsible for the arising of specific correlations among 
the carriers. We remind that determining if a given quantum 
transformation is compatible with a Lindblad structure is in 
general a computationally hard problem [jl Also we notice 
that a Lindbladian structure for the global system S in gen- 
eral may introduce non-Markovian elements in the dynamics 
of the subsystems that compose it, which also are far from 
trivial to characterize, e.g. see Refs. ||12J . To bypass such 
difficulties in our analysis we will thus adopt a rather prag- 
matic approach, deriving the dynamical evolution of S from 
a collisional model llT3lfT4ll in which dissipative effects orig- 
inate from a sequence of weak but frequent interactions with 
a collection of uncorrected particles which mimic the system 
environment. Consider hence a multipartite quantum system 
S, consisting of M - not necessarily identical - ordered sub- 
systems S*!, 5*2, • • • , Sm- In what follows each subsystem is 
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FIG. 1 : Schematic of the process. The horizontal lines describe an or- 
dered set of carriers Si, 82, - ' ' which interact with an ordered set of 
(possibly infinite) identical local sub-environments E\,E2, - ■ ■ via 
local unitaries Us^En- Between collisions each sub-environment 
evolves according to a map M. The overall dynamics can be de- 
scribed as a ordered sequence of row or of column super-operators 
(visualized by the rectangular sets in the figure). 



supposed to interact with a multipartite environment E con- 
sisting of a large number of sub-environments £'2, • • • via 
an ordered sequence of pairwise interactions (for a pictorial 
representation see Fig. [T]). As in |fT3] 114 1 the pairwise colli- 
sion between the subsystem Sm and the sub-envirinment _E„ 
is described by a local unitary Us^e„ = cxp[— zgiJ^^^^^ Ai] 
characterized by a collision time A< and by the intensity pa- 
rameter g, and generated by the Hamiltonian coupling Hs^e„ 
which (without loss of generality) we write as 

Y: Ail « 4! , (1) 

e 

with Ag Hermitian. Accordingly the m-th car- 

rier interacts with the first n elements of the environment £ 
through the joint unitary 

f^slf - Us^.E^ Us,^,E,,_, ■ ■ ■ Us^,E, Us^.E, , (2) 
(the presence of a local free Hamiltonian evolution operat- 
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ing between the collisions can be included in the model by 
passing into the interaction picture representation and replac- 
ing A^g ^ with the corresponding evolved operators). Finally 
to account for the internal dynamics of the environment, we 
assume that between two consecutive collisions each sub- 
environment evolves according to a completely, positive, trace 
preserving (CPT) map A4. 

Consider then the case where the Af subsystems of S are 
initially in a (possibly correlated) state p(0) while the sub- 
environments of £ are all prepared into the same input state 
(which, as in Ref. |13|, represents some equilibrium state of 
the particles of the reservoir). For the sake of simplicity in the 
following we will work under the hypothesis that 
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where we use the symbol {■ ■ ■)x to represent the trace over 
the system X, and 7W" to represent the channel obtained by 
applying m times the map A4. The assumption ^ allows us 
to rigorously define the continuous limit of the model. It is 
worth noticing however that it does not imply any loss of gen- 
erality as it can always be enforced by moving into an inter- 
action representation with respect to a rescaled local Hamilto- 
nian for the system S. 

After the interactions with the first n element of £ the global 
state R{n) of the system and of the environment is obtained 
from the initial state p(0) r;*"" as R{n) = >V("'*^)(p(0) (g) 
7?®"), where yv("'*^) is the super-operator which describes 
the collisions and the free evolutions of £. As schematically 
shown in Fig. [T] it can be expressed as a composition of row 
super-operators stack in series one on top of the other 



(4) 



where Tls,„.£ := ^A'^" o Ug^^ g. Here given, a unitary trans- 
formation U, we define ■ ■) = U{ - ■ ■)W. Also we use the 
symbol "o" to represent the composition of super-operators 
and A^**" to represent Mei ° • ■ ■ ° A^b^, A^b^ being the 
map A4 operating on the j-th element Ej of £. The trans- 
formation Tls^,£ describes the evolution of Sm in its inter- 
action with £ plus the subsequent free evolution of the latter 
induced by the maps A4. Alternatively, exploiting the fact 
that for m' m, n' n the operators Us^,e„ and Us^, ,e^, 
commute, can also be expressed in terms of column 

super-operators concatenated in series as follows: 



yy(«,M) ^ Q 



S,E, 



where for all j = 1 , 



,n. 



S,Ei 



(5) 



Cs,Ej Mej o i^SM,E, o • • • o Mej o ^Si,Ej ■ (6) 

Thanks to Eq. (jSj) we can now write the following recursive 
expression for R{n): 



(7) 



I. THE MASTER EQUATION 

For a particular class of interaction unitaries, the Authors 
of |14| have shown that the collision model leads to a dynam- 
ics which can be described by a Lindblad super-operator via 
direct integration of the equation of motion. Here we intro- 
duce an alternative approach which allows one to derive a ME 
for the reduced dynamics of the many-body system S in our 
generalized multipartite collision model. The details of the 
derivation can be found in the Appendix. We simply assume 
a weak coupling regime where we take a proper expansion 
with respect to the parameters g and At which quantifies the 
intensity and the duration of the single events. In particular 
we work in the regime in which oAt is a small quantity and 
expand the dynamical equation (mi up to 0((gAt)^), i.e. 

Rin+l) = [ls,E„^,+C's^E„+,9At (8) 
+ ClE„^AgAtf]{R{n)^,j) + 0{{gAtf) , 



where Is,e„+i is the identity superoperator while C'g ^ 
and Cg ^ are the first and second expansion terms in gAt 
of the superoperator Cs.e,^+i, respectively. Tracing over the 
degree of freedom of the environment the resulting equa- 
tion defines the incremental evolution of the density matrix 
p{n) :— {R{n))£ of S when passing from the n-th to the 
(n + l)-th collision. The continuos limit is finally taken by 
sending At to zero while g and n explode in such a way that 
nAt and g^At remains finite, i.e. 



lim n At 

At^O+ 



t < oo 



lim g^At 



7 < OO . (9) 



Notice that while the first condition is necessary to properly 
define the axis of time, the second is needed to guarantee that 
S fills the interactions with £. Indeed one easily verifies that 
the linear terms in g do not enter in the dynamical evolution of 
S since (C^^^ ^ (Rin) iSi rj)) £ = due to the assumption (3 1. 

Defining hence p{t) — \im^t_yQ+ p{n) the reduced density 
matrix of S at time t, and f){t) :— lim^t_j.o+ 
time derivative, from Eq. &\ we get the following ME: 



its 



M 



(10) 
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This is mathematically equivalent to the standard derivation 
of a Markovian ME for a system inetracting with a large 
environment, in which one assumes that the overall system- 
environment density operator at any given time t of the evolu- 
tion factorizes as in p{t) ® rj where 77 is the environment den- 
sity operator The two scenarios are however different. In the 
standard case the reason for which the environment state is un- 
changed is because it is big. In our scenario, consistently with 
the collision model, the environment state is constant because, 
as we said, each subsystem collides briefly with a sequence of 
sub-environments all initially in the same state. Of course one 
expects a strongly non markovian behavior if a given subsys- 
tem interacts repeatedly with the same sub-environment |15|. 

The ME ( [TO] i contains both local Lindblad terms (i.e. Lind- 
blad terms which act locally on the m-th carrier) and two-body 
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non local terms which couple the M carrier with the m' > m. 
More precisely the m-th local term is the super-operator 



A. Correlations 
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Equation (13i obeys to proper time-ordering rules which 



where the non negative matrix 7; 



is given by 



7; 



(11) 



(12) 



guarantee that the dynamical evolution of Sm is not influenced 
by the subsystems that follow it in the sequence, while it might 
depend in a non trivial way on the carriers that precede it. 
Indeed when traced over the degree of freedom of the second 



carrier Sm', the cross term T)^^^-^, nullifies, i.e. 



rn'm' V 



= 0, 



(17) 



with 7 as in Eq. ^. Equation ^ defines the correlation ma- ^hile in general it does not disappear when tracing over S, 



(i) (£') 

trix of the sub-environment operators and evaluated 
(for the infinitesimal time interval At) on the density matrix 
which describes the state of the sub-environment 
after to — 1 free evolution steps 1 16|. For m' > m the cross 



terms of Eq. ( 10 1 are defined instead as 



V, 
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with [• • • ,•••]- being the commutation matrix and ^mm' 
ing the complex matrix ifTTI 



(14) 



The coefficients 7^f„'l introduce cross correlation among the 
carriers and depend upon their distance m' — m. Furthermore, 
similarly to the the terms of Eq. ( [T2) l, they also depend on 
TO — 1 due to the fact that the model admits a first carrier 
However if we assume that for large to the sequence A4"^{ri) 
converges to a final point rjQ, then we can reach a stationary 

(i £') 

configuration where (for to ^ 1) 7„j'„// only depends upon 
the distance m' — to while 7™'^ becomes constant in m, i.e. 



7m, m' 

I m 



(15) 
(16) 



A similar behavior is obtained also if we assume 77 to be a 
fix point for Ai (a reasonable hypothesis if £ is supposed to 
describe an environment in its stationary configuration). In 
this case Eqs. (15 1, (16 1 hold exactly for all to and m', with 770 
being replaced by 77. Finally a case of particular interest is the 
one in which Ai is the channel which sends every input state 
into 7] (this is the extremal version of the last two examples). 
Under this condition one expects that no correlations between 
the various carriers can be established as the environmental 
sub-systems are immediately reset to their initial state after 
each collision. Indeed in this case we have A4{9) — {0)e V 

for all operators 6, which, thanks to Eq. I3 1, yields J^'m' = 

j{Bpr,)E{B^E''v)E = and hence P^J - 0. 



(it does disappear however if all the coefficients 7^'^/ are 



real, see below). The evolution described by Eq. ( lOi is thus 
non-anticipatory fTSl, or in the jargon introduced in Ref. [19|, 
semicausal with respect to the ordering of the channels uses. 
To see this explicitly consider the evolution of the reduced 
density matrix pi.2{i) of the first two carriers obtained by tak- 
ing the partial trace of Eq. ( [T0| over all elements of S but 
Si and Noticing that ' ' — ™d exploiting 



Eq. ( 17 1 we get 



pi,2(t) = A(pi,2(t)) + /:2(pi,2(t)) +2?i;^^(pi,2W) -(18) 

The resulting dynamics is purely Markovian in full agree- 
ment with the fact that 5*1 , S2 couple weakly and sequen- 
tially with sub-environments £ which have not interacted yet 
with other carriers. Tracing over S2 we can then derive the 
dynamical equation for Si, i.e. pi{t) — Ci{pi{t)), which 
again is Markovian. Vice-versa the dynamics of S2 cannot be 
expressed in terms of a close differential equation for P2{t) 
alone. Indeed by taking the partial trace of Eq. ( 18 1 over Si 
we get 



^2(p2(i)) 

27^Im[7(' 



(19) 



(«')i 
2 J 



4„(A5,W,Pl,2(t))i 



where the last term expUcitly depends upon the joint density 
matrix of S2 and Si fSOl . This formally shows that in general 
Si acts as controller for S2, while no back-action is allowed 
in the model. 

A case of special interest is represented by those situa- 

(i i') 

tions in which the matrices 7„'„/ are real. When this hap- 
pens also the partial trace over Sm of 'D^mm' nullifies, i.e. 

(^mm' (■ ■ ' = 0- Accordingly the evolution of any sub- 
set of S is independent from the evolution of the remain- 
ing carriers. In this case hence our model becomes non- 
anticipatory with respect to all possible ordering of the carri- 
ers, describing hence a non-signaling evolution 1 19] in which 
the reduced density matrix of each carrier evolves indepen- 
dently from the others. For instance in Eq. ([19]) the second 
line disappears yielding a Markovian equation also for P2{t), 

i.e. p2{t) = L2{p2{t)). 

a. Example:- As an application we focus on the case in 
which the carriers and £ form two sets of independent bosonic 
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modes. In particular defining a,„ and 6„ to be annihilation op- 
erators of the modes Sm and En respectively, we consider the 
Hamiltonians Hs^^^e„ — o,m^bl^+al^(Sibn. We also take 77 as 
the vacuum state of En- and as a lossy Bosonic quantum 
channel of transmissivity k. Notice that with these choices the 
Hermitian operators and entering in Eq. |l| are just 
quadrature operators of the fields, and that Eq. (j3]l is automati- 
cally verified for all m since J\4 (77) — rj. The resulting model 
describes a correlated quantum channel analogous to that of 
Ref |8| which mimics the transmission of a sequence of op- 
tical pulses along an attenuating optical fiber characterized 
by finite relaxation times. The corresponding local • • ) 
and cross term 2?^^, entering in the final ME |lo| become 
respectively | {2a„(- • • )a„i - al^a^i- ••)-(••• )a^„a„i} 

-[(• ••)aT^,a,n']-} which ex- 



I J. 

and7K 2 {[a„(. . . ),aT„,]_-[^- • • ja;„,a„ 

hibit an attenuation of the signals and an exponential decaying 

in the correlations (in particular ' ' ) coincides with 

the cross term derived in Ref. |21 1 for a collection of QED 

cavity modes coupled in cascade). 



II. CONCLUSIONS AND PERSPECTIVES 

In deriving the ME ( [T0| we assumed a specific order- 
ing for the carriers of the model which implies that each 
elements in the sequence 5i,S'2,--- ,Sm can influence 
only the dynamical evolution of those which follow. This 
assumption was specifically introduced to account for the 
causal correlations that are present in many memory quantum 
channel models |18|. The collisional model however can be 
generalized to include more general correlations. For instance 
cyclical correlations can be accounted by identifying Si with 
the {M + l)-th element of the set of carriers in such a way that 
Sm can influence its dynamics. To do so it is sufficient to add 
an independent set F of sub-environments Fi, F2, • • • , F^v 
which couple with iS following a new ordering in which (say) 
all the carriers are shifted by one position (i.e. the element of 
F first interact with 5*2, then with 6*3, 5*4, • • • , Sn, and finally 
with ^i). A part from the new ordering the new couplings are 
assumed to share the same properties of those that apply to 
£ (in particular we require that identities analogous to those 
in Eqs. Q, (|9]l hold). Under these conditions (and assuming 
no direct interaction between £ and F) the ME ([TO]l will 
acquire new extra terms which directly couple each carrier 
with all the others. Specifically given m' > m we will 
have both a standard contribution of the form 2?^'^^^, as in 



Eq. (10 1 but also a contribution in which the role of m and 
m! are exchanged (i.e. something like 2?^^^^^) that originates 
from the couplings with F. From this example it should be 
clear that by increasing the number sub-environmental sets 
and by properly tuning their interactions with S any sort of 
correlations can be built in dynamical evolution of the system. 

Appendix A: Technical sections 



In this section we give the detailed derivation of Eq. (10 1 
and discuss its generalization to the case of non uniform col- 



lisional events. Subsequently we show how to include free 
evolution terms induced by local Hamiltonians operating on 
the carriers in the derivation of the ME. 



1. Derivation of Eg. (lO) 

The starting point of the derivation is Eq. (j9]) which under 
partial trace over £ yields the identity 

p{n + 1) - p{n) + igAt) (c^,b,.+i (^(^) ® v)) ^ (Al) 
+ (.gAt)2 {c'Ie„+^ (Rin) ® ry))^ +0{{gAtf) , 

In this expression we need to specify the super-operators 
and obtained by expanding Cs.e^+i up to 

the second order in gAt. To do so we notice that for each m 
and j, the super-operators Us^,Ej admit the following expan- 
sion, 

+0{{gAtf) , (A2) 
with Is„^,Ej being the identity map and with 



^S^.Eji---) ■= Hs,„,E,{- ■ ■)Hs^^Ej 
^^Sr^.EjA- ■ ■ ) 



(A3) 



(A4) 



where [• • • ,•••]- and [• • • ,•••]+ represent the commutator 
and the anti-commutator brackets respectively. From Eq. (|6]l 
it then follows that 



M 



^S,Ej 



r" 



Y.^Er^'-^L.E,oM'S-\ (A5) 

(A6) 



771 — 1 



p it-'. 



with 



■"S.Ej 



EkyiM — m+l Till „ A,*"* — 1 



m—1 

M M-1 



m'— m+l rn— 1 



(Al) 



Replacing this into Eq. (jATJ we first notice that due to Eq. ([3]) 
the linear term in gAt nullifies. Indeed we get 
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m I 



^EE[41'^(") 



m £ 



(A8) 



Vice-versa for the second order terms in gAt we get two contributions. The first is 



- I E E(«^ [242p(")4i - 4142p(j^ - 



(A9) 



with C„i as in Eq. ( [TT) . The second term instead is 



m' >m 



A' 



(^) 



7 , 

m' >m 



with as in 



Eq. |l3| . Replacing all this into Eq. ( |ai[ i gi' 



p(n + 1) - _ g^At 



At 
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m' >m 



(AlO) 



(All) 



which enforcing the limit (|9]l yields the ME (10 1. 

It is worth noticing that the above derivation still applies 
also if the collisional Hamiltonians ([TJ are not uniform. For 
instance suppose we have 



H 



(n,l) 



B 



(A12) 



where now the operators acting on the carrier Sm are allowed 
to explicitly depends upon the n index which label the colli- 
sional events, and similarly the operators acting on the sub- 
enviroment are allowed to explicitly depends upon the index 
771 which labels the carriers. Under these conditions one can 



verify that Eq. (Alli still apply. In this case however, to 
account for the non uniformity of the couplings, the condi- 



tion Q needs to be generalized as follows 



(77)) =0, V777,^ 



Furthermore both £„, and entering in Eq. (All 



E 



(A13) 



become explicit functions of the carriers labels ana of 
the index 77, which plays the role of a temporal parame- 
ter for the reduced density matrix ^(77). Specifically the 
new supe r-ope rators are still defined respectively as in 



Ag-* and with the coefficients {B'^^'B^ 

M 



Eqs. (A9l and (AlOi with the operators ^^"^^'^^ instead of 



^E 



'{B 



placed by {B 



{B\ 



E 
\^E 



>B 



M"'-\r,))E 
M"'~Hv)) )e re- 
E M"'-\f])}E and 

A^™^^(77)) ) E respectively. 
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The continuos limit (19]) can also still be defined by identi- where iJ^"^ ^ is the unitary that defines the colHsions of S,, 



fying limAt^o+ A 



(n- 



with the element A^P 



(t) of a one 

parameter family of operators. As a result we get a time- 
dependent ME characterized by a Lindblad generator which 
explicitly depends on t. 



2. Including local free evolution terms for the carriers 

Assume that between two consecutive collisions, the carri- 
ers undergo to a free-evolution described by a (possibly time- 
depedent) Hamiltonian Hs{t) :— J2m ^Smi^) which are lo- 
cal (i.e. no direct interactions between the carriers is allowed). 



Under these conditions Eq. ( 10 1 still holds in the proper inter- 
action picture representation at the price of allowing the gen- 
erators of the ME to be explicitly time dependent. 

To see this we first notice that under the assumption that 
the collision time At is much shorter than the time interval 
that elapses between two consecutive collisional events (i.e. 
At <C T„ — T„_i), the unitary operator which describes the 
evolution of the m-th carrier in its interaction with £ is now 
given by 



(n) 



) 



s,„,£;„_i 



(A14) 



where Us„^.e„ are the collisional transformations, r„ is the 
time at which the n-th collision takes place, and where 
Vs^{Tn,Tn^i) ■= Texp[-i J^^\^ dt' hs^{t')] is the uni- 
tary operator which describes the free-evolution of Sm be- 
tween the {n — l)-th and the n-th collision (in this expression 
Texp[- • • ] indicates the time-ordered exponential which we 
insert to explicitly account for possibility that the hg^^^ will be 
time-dependent). Define hence the operators 



A 



(A15) 



and the Hamiltonian 



(n.i) 



B 



(A16) 



with the sub-environments in the interaction representation, 
i.e. 



U. 



with 



Srr^■£ 



Us 



(A18) 



(A19) 



Similarly we can express the super-operators )/v("'^^) as 

y^(n,M) ^ V5(r„,0)o yy("^^'^) , (A20) 
y^KM) _ Cs,E„o---oCs,E^ , (A21) 

Cs.E, ■■= Me, o Us^j,e, o • • ■ o Me, o Us,.Eih'2-2) 

with V5(t„,0) being the super-operator associated with 
the joint free unitary evolution obtained by combining 
all the local terms of the carriers, i.e. V5(t„,0) := 
Vsi {Tm^) - ■ ■ Vsm {'Tn, 0). Defining hence R{n) the state of 
S and of the first elements of £ after n collisions in the inter- 
action representation induced by Vs{Tn, 0) as 



= yj(T„,0) i?(n) Ks(t„,0) , 



(A23) 



we get a recursive expression analogous to Eq. (j7]l with 
Cs,e„+i replaced by C5.£;„+i, i-e. 



R{n + 1) = Cs,E„+,{R{n)(g)ri) 



(A24) 



which describes the coupling between S„i and £ in the inter- 
action representation associated with the free evolution of Sm. 

Notice that the operators A^!^^^ are explicit functions of the p(t) obeys^to a ME analogous to^Eq^lfTol with the opera 



More precisely this expression formally coincides with that 
which, as in the case described at the end of the previous 
section, one would have obtained starting from a collisional 
model in which no free evolution of the carriers is allowed but 
the collisional events are not uniform. Indeed the generators 
of the dynamics Hs,^.e„ do have the same form of the 
Hamiltonians ( |A12| l. Following the same prescription given 
there, we can then get an expression for the reduced density 
matrix p{n) ~ {R{n))£ which represents the state of the 
carriers after n collisions in the interaction picture with 
respect to the free evolution generated by Hs{i). Enforcing 
the limit (j9]) under the condition (A13i, one can verify that 



index n which labels the collisions as in the case of Eq. ( A12 1 



(here however the terms operating on £ are kept uniform). Ob- 
serving that for all t one has Vg^ (r^, Tj-ijV s^ {ti~i, ti~2) 
{Tf., Tg^2) we can now write Eq. ( A14 1 as 



tors Ag ' being replaced by the time-dependent operators 



► 0+ 



r(n/) 



(«) . 



Vs„(t-„,0) U, 
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